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Abstract. We analyze the d'Alembert equation in the Godel-type spacetimes 
with spherical and Lobachevsky sections (with sufficiently rapid rotation). By 
separating the t and £3 dependence we reduce the problem to a group-theoretical 
one. In the spherical case solutions have discrete frequencies, and involve spin- 
weighted spherical harmonics. In the Lobachevsky case we give simple formulas 
i for obtaining all the solutions belonging to the sectors of the irreducible 

unitary representations of the reduced Lorentz group. The wave equation enforces 
restrictions on A and the allowed (here: continuous) spectrum of frequencies. 
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1. Introduction 



When a rapidly rotating gas of bosons is cooled down below the Bose-Einstein 
condensation temperature a regular lattice of vortices (Abrikosov lattice) emerges. 
The velocity vector field which initially circulated around the only axis of rotation 
now rotates locally around every vortex core. Interestingly, in the General Relativity 
Theory there also exist both types of rotation. A homogeneous spacetime with the 
latter type of rotation was discovered by Godel pQ long time ago. Chapline and Mazur 
have argued recently [2] that a BEC with an array of vortices and a Godel-type metric 



will provide a solution of the Einstein-Gross-Pitaevskii system (in the coarse-grained, 
v <C c s limit). By Godel-type we mean homogeneous stationary spacetimes in which 
matter rotates around every point relative to the "compass of inertia" [21 H| . As in 
isotropic cosmology, the hypersurfaces H =const" can have either Lobachevsky, flat 
(3JT}| or spherical geometrjQ with the radius of curvature being a parameter. The angular 

velocity of rotation provides the second parameter of the Godel-type spacetimes. 

In this paper we set ourselves a modest task of analyzing the d Alembert equation 
in the Godel-type spacetimes with spherical and Lobachevsky sections (in the latter 
case we confine ourselves to the case of sufficiently rapid rotation). We shall, of course, 
profit from the work of a number of authors who have, from different angles, attacked 
this problem before, El El El E3] ■ Therefore, our goal here is not to discover, but 
rather to consolidate, simplify (and where necessary - to correct) the work of others 
in order to lay solid foundations for further progress. 

This short paper has one main section (the second one), which is divided into 
two parts treating the dAlembert equation in the Lobachevsky and spherical cases 
separately. We have also included appendices on the Lobachevsky geometry and on 
the sources of Godel-type spacetimes for quick reference. 

I The type of geometry is controlled by, among other things, the sign of the cosmological constant. 
The spherical one corresponds to a positive vacuum energy currently supported on the cosmological 
scale by multiple observations. 
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2. Solving the wave equation by an Ansatz 

In order to find the solutions of the d'Alembert equation in the spacetimes of Godel 
type we will make use of the fact, that these spacetimes are stationary homogeneous 
and axially symmetric. Thus, we will assume that the solutions, 'J, of the wave 
equation, = 0, can be separated according to 

*(t, x, <p, x 3 ) = e - iut e ik ^ ■ ip(x, ip), (1) 

where the function tp(x,tp) should be periodic in p with the period of 27r. In the 
following we shall classify the allowed functions, ip(x,p), in the Lobachevsky and 
spherical cases separately 



2.1. Lobachevsky case 

The original Godel spacetime is a stationary four-dimensional spacetime with the 
t =const hypersurfaces being just H2 x R, where H2 is the Lobachevsky space (or 
hyperbolic plane, see appendix A). These hypersurfaces are threaded by a constant, 
homogeneous gravimagnetic field perpendicular to the Lobachevsky space (i.e. along 
the X3 direction). There family of Godel- type metrics pQ is parameterized by the 
"amount of rotation" £1 and the magnitude of Lobachevsky curvature denoted by I 
(see appendix B): 

d(2ls) 2 = (dt + Hdp) 2 - dx 2 - D 2 dp> 2 - dxl, (2) 

where all the coordinates are dimensionless, and 

H = 2asinh 2 (x/2) (3) 

L»=sinh(a;), (4) 

with a = £1/1. Both t and £3 range over the whole R, while x £ [0, 00) and 
(p £ [0, 2tt) with the usual identification of and 2tt. The metric is of the form 
dr 2 — (dt+Aidx % ) 2 — hijdx t dx\ with the Lobachevsky space metric hy and a constant 
gravimagnetic (Lense-Thirring) field on this space 

B> = = 2a, (5) 

with the same constant being equal to the contravariant X3 component of the rotation 
vector of the (d t ) a vector field (velocity field of the sources). The case O = 
corresponds to a static spacetime with Lobachevsky spatial sections, while il 2 = 2l 2 
corresponds to the original Godel metric. 

In this paper we deal only with spacetimes of sufficiently rapid rotation £1 > I 
(that is, a > I). In order to solve the wave equation in these spacetimes, we will make 
use of the underlying Lobachevsky structure. The d'Alembert equation for \& is 

□* = [1 - a 2 tanh 2 (x/2)] *« + ° — tt ty - 

cosh (x/2) 

- . \. r ^yy - ^33 ~ ®xx ~ COth^)*^ = 0. (6) 

smh (x) 

The following three operators, which are Killing vectors of the Godel spacetimes 
multiplied by —i, fulfill the same commutation relations as the Killing vectors Li 
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of the Lobachevsky plane (see appendix A): 

J\ = —ismtpdx — i coth(a;) cos^)^ — iacosipt&nh(x/2)dt, (7) 

J-2 = i cos ipd x — i coth(x) sm((p)d v — ia sin ip tanh(a;/2)<9 t , (8) 

Jo = - id v - iad t , (9) 

while the remaining two Killing vectors, dt and 83, obviously commute with the above 
ones. We notice further, that making use of the separation ansatz ([l} for 'J, we are 
lead to 

□ * = -C 2 * - (1 -a 2 )uoH- ^33, (10) 

where here: C2 = Jq — Jf — Jf ■ The wave equation therefore reduces remarkably to 

C 2 ^{x, <p) = [{a 2 - l)uo 2 + k 2 } ip) (11) 

Because we are interested in the regime a 2 > 1, it remains to determine the possible 
eigenfunctions, ip(x,ip), with a positive eigenvalue of the Casimir operator C2 which 
was defined above. Because of the ansatz, the operator C2 becomes now a second 
order partial differential operator acting on ?/>(a;, ip). 

In order to determine the eigenfunctions, let us now regard C2 (and, in fact, the 
whole set J a ) as an operator acting in the Hilbert space of square- integrable functions 
of x, tp with the scalar product provided by the Lobachevsky product, (|A.4|) . 

It is helpful to recall the classification of the unitary irreducible representations of 
the reduced Lorentz group £3 (two boosts and one rotation) and of its corresponding 
spinor group 63 =SU(1,1) of Bargmann lllj. For positive eigenvalues of C2 one is lead 
to what Bargmann calls discrete classed, and D^. Within the eigenvalue of 
C 2 is A(A =F 1) and the eigenvalues of Jo start with A and increase (without bound) 
by one with each application of J+ = Ji + iJ 2 in the class (or, respectively, start 
with — A and decrease (without bound) by one with each application of J_, for the 
class) . In the number A is real and positivcfjj], while in it is real and negative. 

Let us determine all these functions explicitly. For the complex parameter 

z = tanh(.T/2)e l¥ \ (12) 

we note that 

d v = izd — tzd, 
where d = c\. We find 

J+ = z 2 d-d + az, (13) 
J_ = d - z 2 d + az, (14) 
Jo = zd — ~zd + a, (15) 

where 

a = —auj (16) 
and J± = J\ ± iJ%. We may verify, that 

[J+,J-}= -2J , (17) 
C 2 = J (Jo + 1) - J-J+, (18) 
C 2 = Jo(Jo-l)- J+J-, (19) 

§ Note, that our convention on the sign of the Casimir operator, C2, is different from that of 
|| Bargmann notes (see §4 of that the universal covering group £ covers S3 infinitely often. He 

restricts himself to single- valued representations of S3, and as a consequence treats only integer and 
half-integer A's. We see no reason for such a restriction here, and therefore leave A real. 
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which are all standard relations in the theory of representations of SU(1,1). 

For a given A the families of simultaneous eigenfunctions of C2 and Jo are 
constructed upon a reference vectors, denoted here by ip^ (or V'o"); annihilated bjfij 
J_ (or J + ). This leads to two families, and D7. 

2.1.1. Structure of . We find immediately 

= (1 - zz) x z a -\ (20) 

where A is the eigenvalue of Jo, and A(A — 1) the eigenvalue of C*2- The function 
ipQ depends, via the last term, on the angle cp. Since it must be periodic in ip, with 
identical values for tp = and (p — 2tt, the number a — A = m must be integer. In the 
original variables we have 

^+ = cosh(a;/2)- 2A - m sinh(a;/2) m e' lmip . (21) 

Further restrictions, 

m>0, A>i, (22) 

are obtained from the requirement of square integrability w.r.t. the product (|A.4[) at 
x = and at infinity, respectively. 

A connection with the methods of ordinary differential equations is provided by 
the following consideration: the eigenvalue problem for the Casimir operator, 

C 2 ^j = A(A - 1)V, (23) 

with A > 1 leads, with a separation ansatz ip — f (x) exp(—irrup) , Z 3 m > to the 
equation 

2a(a + m) m 2 



w(w-2)f"+2(w-l)f+ 



/ = A(A-1)/, (24) 



w(w 

where w = 2 cosh 2 (x/2). This is an equation with three regular singular points, at 
(0, 2, 00), with the exponents (a + a — (^, — y) and (A — 1, -A) respectively. 
Evidently, the function(s) we have found by algebraic methods, tpQ , have the right 
asymptotics at w = 2 (V'J" ~ {w — 2) "2") and at infinity (i^o ~ w~ x ) to be square 
integrable at these extremes. 

Let us now summarize the results. The solutions are characterized by the 
momentum in the X3 direction (^3) and by the representation parameters (A, k) with 
k G {A, A + 1, . . . , 00}. The wave equation has the form 

A(A - 1) = (3 2 a 2 + k 2 , (25) 

where f3 2 = 1 — 1/a 2 > depends only on the parameter a of the spacetime, and 
a — A > is an integer. The above dispersion relation is unusual by the sign in front 
of k 2 . The solution for a is given by 
1 



-^(A-l)-^ 2 , (26) 

but only if a > A, which is possible only if A > a 2 . Consequently, because of (TlGl) . 
there also exists a gap in frequencies: \lo\ > a. 

% The necessity for such a vector is clear; as in the case of the rotation group - its absence would 
lead to negative- norm states created by successive application of either J_ (or J+). 
+ Even though on the surfaces of constant t and £3 the dip vector becomes time-like for large x, 
allowing for functions growing exponentially with x (omitting the second restriction) would certainly 
make no sense. 
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2.1.2. Structure of D x . The function 

% = {l-zz)- x z x - a , (27) 

is an eigenfunction of C2 to the eigenvalue A(A + 1), and of Jo to the eigenvalue A. It is 
square-integrable if A < —1/2 and if m = A — a (which is an integer) is non- negative. 
Consequently, C2 is non- negative and ipQ is square-integrable if a < A < — 1. Because 
of that, the frequency ui is now positive, and the wave equation 

A(A + 1) = (3 2 a 2 + k 2 , (28) 

has solutions 

a = --^\{\ + l)-kl (29) 
only if A < —a 2 (in order to satisfy a < A). Here, again \ui\ > a. 

2.2. Spherical case 

The case of spacetimes with spherical (x, if) hyperplanes (threaded by a homogeneous 
gravimagnetic field) is simpler, and the resulting functions, ip(x, if) belong to a known 
family of functions (spin "a" spherical harmonics). In this case, H = 2a sin 2 (a;/2) 
and D = sin(x) (where x S [0, 71") has the interpretation of the latitude). The wave 
equation can be solved similarly, though in this case the d'Alembert operator is given 

by 

□ * = [C 2 - (1 + a 2 )u 2 + fcf] * = 0, (30) 

where C2 is the Casimir operator of a modification of the algebra the SO(3) of 
(symmetry) generators: 

L+ = z 2 d + d + a Zl (31) 
L_ = -d - z 2 d + az, (32) 
L3 = zd — zd + a, (33) 

and 

z = tan {x/2)e lip . 
The above operators fulfill the SO (3) relations: 

[L+,L_]=2L 3 , [L ± ,L 3 ] = T L±. 

The representation spaces are grouped in sectors, each with a definite value of 
L 2 = A(A + 1) and L 3 e [—A, A], where A is a positive integer or half-integeiQ. 

The lowest states, ip-\, which are annihilated by L_ with the eigenvalue —A of 
L3, and the highest states ipx (annihilated by L + with the eigenvalue A of L3) can 
now be determined. One finds 

4'-x = (l + zz)- x z x+a , (34) 
iPx = {1 + zz)- x z x - a . (35) 

* Only the integer or half integer values of A are allowed, for otherwise one would be lead to standard 
contradictions (eg. states of negative norm, after sufficiently many applications of L+). 
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These function must be periodic in ip, with identical values for ip — and ip = 2tt, and 
therefore a + A = m must also be an integer. Now it is easy to see, that either there 
holds 

m > 0, (36) 
m < 2A, (37) 

or the function ip-\ will not be square integrable on the sphere (in x, tp) because 
of singularities at x — or x = ir respectively. Thus, for a given value of 
A G {i, 1, . . . , oo} there remains a freedom to choose m € [0, 2A]. This can be viewed 
as a condition specifying a, so that the structure is the following: 

There are infinitely many sectors with integer or half-integer a. In each of these 
sectors there are eigenspaces of L 2 to the eigenvalues A(A + 1) with A > \a\ (with an 
integer A + a). Each of such eigenspaces is spanned by 2 A + 1 functions generated by 
a successive application of L + to ip—\. 

The eigenspace A = i, a = — ^ consists of two functions first encountered by 
Dirac [12] in the context of magnetic monopoles, 

V^o = cos(x/2) for k = ~, (38) 

ip x = - sm{x/2)e iip for k = -. (39) 

The functions we have found are also known as the monopole harmonics or the spin-a 
spherical harmonics 114], so that further properties (eg. completeness for each a) 
and identities among these functions are known. 

We may now summarize the results taking into account the wave equation. Each 
solution is characterized by the parameters: frequency (u>), momentum in the X3 
direction (ks), sector number a (which is an integer or half- integer), and the pair 
(A,fc) with A > |o| and k € [-A,A]. With 

uj = - (40) 
a 

the frequency is fixed by the sector number a. Such a quantization of frequencies has 
an alternative, kinematical (gauge) derivation pointed out long time ago by Mazur 
[15] . The only role of the wave equation is to provide a relation between £3, A and a: 

A(A + l) + fe| = (l + l/a 2 )a 2 . (41) 

Thus, in each sector (for each energy) fc 3 belongs to a discrete bounded set. 



3. Conclusions 



Starting with the separation ansatz ([T]) we have solved the d'Alembert equation in 
the spacetimes of Godel type. In the spherical case the functions with the frequency 
uj = n/2a, n G Z turned out to coincide with the n/2-weightcd spherical harmonics. 
In the Lobachevsky case we have given simple formulas for obtaining all the solutions 
belonging to the sectors of irreducible unitary representations of the reduced 
Lorentz group. The wave equation enforced a 2 < |A| e R and a gap in frequencies 
uj G R\(-a,a). 
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Appendix A. Preliminaries on the Lobachevsky space 

The 2-dimensional Lobachevsky space (known also as the hyperbolic plane, H 2 ) is a 
space defined as the hypersurface t 2 ~x 2 — y 2 of the 1+2 dimensional Minkowski space. 
The intrinsic metric of this space is given by 
ds 2 = dr 2 + sinh 2 (r) dtp 2 , 

with r G [0, oo), f € [0, 2n). (This map covers the whole space.) The three Killing 
vector fields, 

Ki = sm(ip)d r + coth(i/>) cos(<p)d v (A.l) 

K 2 = — cos((f)d r + coth(f/') sm(ip)d v (A. 2) 

K = d v , (A.3) 

by complete analogy with rotations on a sphere, when multiplied by — i lead to the 
generators of the respective symmetries on the Hilbcrt space of square integrable 
functions on H 2 , with the standard product 

(/, g) = j Jg sinh(r) dr dip. (A.4) 

These generators, denoted here by L,, fulfill the commutation relations of the algebra 
of SU(1,1) generators: 

[L ,ii] =*i 2 (A.5) 

[L 2 ,L ]=zLi (A.6) 

[L 1 ,L 2 \ = -iL , (A.7) 
and commute with the quadratic Casimir operator 

C 2 = L 2 -L 2 - L\. (A.8) 

The group SU(1,1) is not compact, and therefore the representation theory is richer 
than that of the rotation group. Depending on the sign of the eigenvalue of C 2 there 
are, discrete or continuous series of eigenfunctions. Let us finally notice, that 

C 2 = d 2 + coth(rR + , * 3% = ^=di (h^ Vhdj) (A.9) 
smh (r) v/i v J 

is just the Laplace-Beltrami operator on the Lobachevsky space, and therefore here 

C 2 < 0. 
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Appendix B. Sources 

Appendix B.l. Spherical case 

Let us choose the orthonormal tetrad 

e° = dt + mdip, 

e* = Dd<p, 

e R = dr, 

e 3 = dx 3 , 

where the functions m and D are as in section 12.21 The spherical Godel spacetime 
can be generated by [THE] a (neutral) dust, supplemented by a cosmological constant 
(with a positive vacuum energy) and an (source-free) electromagnetic field directed 
along the x% axis: 

F m = Ecos(ax 3 ), 
F 30 = Esm(ax 3 ). 
The Einstein field equations 

G A b = k(pu a u b + T AB ) - Atjab (B.l) 

with k = 8wG/c 4 , and the dust four-velocities, ua, coincident with the (dt) a (and 
with eg) vector fields. We get 

G AB =diag(l + 3A,A,A,-l-A), 

Tab =|Uiag(l, 1,1,-1), 
-Atjab = -Adiag(l, -1,-1,-1), 

where A = a 2 /4. After an additional rescalling of the line element ds 2 ^>R 2 ds 2 , with 
R = 1/21 being the radius of the sphere, we find 

A =-2^< 



P 



16itGR 2 ' 

E 2 = -^(2 + a 2 ). 
AGR 2 ' 

Alternatively, by introducing the vacuum energy, po = —Ac 4 /8irG > 0, we may write 
the last two equations as 

P = a 2 p , 

E 2 = iwp (2 + a 2 ). 

Thus, in the spherical Godel-type spacetime the vacuum energy controls the scale R, 
the dust's density the rotation, while the energy density of the electromagnetic field 
needs to be of the order of the vacuum energy. 
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Appendix B.2. Lobachevsky case 

We use the same tetrad form, but now with the functions m and D from section 12.11 
We get 

Gab = diag(-l + 3A, A,A,1 — A), 

Tab = |Uiag(l, 1,1,-1), 

and the Einstein equations lead to 
1 

~ 2i? 2 ' 



16irGR 2 ' 



For a 2 — 2 we obtain the Godel's original solution, which is generated by the dust 
and the cosmological constant alone. 



